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The gluon propagator in Landau gauge is calculated in quenched QCD on a large 
(32 s X 64) lattice at (3 = 6.0. In order to assess finite volume and finite lattice 
spacing artefacts, we also calculate the propagator on a smaller volume for two 
different values of the lattice spacing. New structure seen in the infrared region 
survives conservative cuts to the lattice data, and serves to exclude a number of 
models that have appeared in the literature. 


1 Introduction 

The infrared behaviour of the gluon propagator is important for an understand¬ 
ing of confinement. Previous conjectures range from a strong divergence Era to 
a propagator that vanishes in the infrared cM 

Lattice QCD should in principle be able to resolve this issue by first- 
principles, model-independent calculations. However, lattice studies have been 
inconclusive up to now0H since they have not been able to access sufficiently 
low momenta. The lower limit of the available momenta on the lattice is given 
by g m in = 27 t/L, where L is the length of the lattice. Here we will report results 
using a lattice with a length of 3.3 fm in the spatial directions and 6.7 fm in 
the time direction. This gives us access to momenta as small as 400 MeV. 


2 Lattice formalism 

The gluon field A M can be extracted from the link variables U^(x) using 
Upix) = e^W/ 2 ) + 0(a 3 ). 

Inverting and Fourier transforming this, we obtain 

A,{q) = Y. e ~ lHX+tl/2) Mx + N 2 ) 


X 

e ~iq^a/2 

2 igoa 




Ar (U^(q)-Ul(-q)) 


( 1 ) 


( 2 ) 
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where £7 M (g) = J2 x e l ^ x U R (x) and A^(q) = t a A a ^(q). The available momen¬ 
tum values q are given by 

q^ = 2irn ll /(aN ll ), n M = 0,..., IV M - 1 (3) 

where N fl is the number of points in the /i direction. The gluon propagator 
D^ b u (q) is defined as 

D^(q)=(A;(q)A b v (-q))/V. (4) 

In the continuum Landau gauge, the propagator has the structure 

D; b (q)=S ab ( 6 ,„-^)D(q 2 ). (5) 

At tree level, D{q 2 ) will have the form 

£><°V) = 1/q 2 . (6) 


On the lattice, this becomes 

£ (0) (g) = l/W^sin^) . (7) 


Since QCD is asymptotically free, we expect that up to logarithmic corrections, 
q 2 D{q 2 ) —> 1 in the ultraviolet. Hence we define the new momentum variable 
Q by 


2 . q^a 
q, n = - sin—, 


( 8 ) 


and work with this throughout. 

The (bare) lattice gluon propagator is related to the renormalised contin¬ 
uum propagator D R (q; /z) via 


D L (qa) = Z 3 (n, a)D R (q; n). 


(9) 


The renormalisation constant Z 3 ( 11 , a) can be found by imposing a momentum 
subtraction renormalisation condition D R (q)\ g 2— M 2 = -K . 

The asymptotic behaviour of the renormalised gluon propagator in the 
continuum is given to one-loop level by 

Dn(q) = D haie {q)/Z 3 = ^{l br(g 2 /A 2 )) ° (10) 

with 

39 —C —41V/ 13 

D 4(33-2 N f ) 44’ 1 ’ 

where both the gauge parameter £ and the number of fermion flavours Nf are 
zero in this calculation. 
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Name 

P 

a* (GeV) 

Volume 

-^conf 

&(?Max 

\9M 

Small 

6.0 

1.885 

16 3 x 48 

125 

2tt/4 

< io- e 

Large 

6.0 

1.885 

32 3 x 64 

75 

2tt/4 

< 10“ 6 

Fine 

6.2 

2.63 

24 3 x 48 

223 

2tt/4 

< 10" 6 


Table 1: Simulation parameters 



Figure 1: Componentwise data, for the small lattice (left) and the large lattice (right). The 
filled triangles denote momenta directed along the time axis, while the filled squares denote 
momenta directed along one of the spatial axes. 


3 Simulation Parameters, Finite Size Effects and Anisotropies 

We have analysed three lattices, with different values for the volume and lattice 
spacing. The details are given in table |l]. In the following, we are particularly 
interested in the deviation of the gluon propagator from the tree level form. 
We will therefore factor out the tree level behaviour and plot q 2 D(q 2 ) rather 
than D(q 2 ) itself. 

Fig. 0 shows the gluon propagator as a function of qa for the small and 
large lattices, with momenta in different directions plotted separately. For 
low momentum values on the small lattice, there are large discrepancies due to 
finite size effects between points representing momenta along the time axis and 
those representing momenta along the spatial axes. These discrepancies are 
absent from the data from the large lattice, indicating that finite size effects 
here are under control. 

However, at higher momenta, there are anisotropies which remain for the 
large lattice data, and which are of approximately the same magnitude for the 
two lattices. In order to eliminate these anisotropies, which arise from finite 
lattice spacing errors, we select momenta lying within a cylinder of radius 
A qa = 2 x 27r/32 along the 4-dimensional diagonalsB 
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4 Scaling behaviour 

Since the renormalised propagator D R (q ; n) is independent of the lattice spac¬ 
ing, we can use ([)]) to derive a simple, q-independent expression for the ratio 
of the unrenormalised lattice gluon propagators at the same value of q: 

Dc(qa c ) = Z 3 (n,a c )D R (q;n)/al = Z^_ ^ 

Df{qa f ) Z 3 (n,af)D R (q- lf j,)/a 2 f Z f a 2 

where the subscript / denotes the finer lattice {(3 = 6.2 in this study) and 
the subscript c denotes the coarser lattice {(3 = 6.0). We can use this relation 
to study directly the scaling properties of the lattice gluon propagator by 
matching the data for the two values of (3. This matching can be performed 
by adjusting the values for the ratios R z = Zf/Z c and R a = a,f/a c until the 
two sets of data lie on the same curve. 

We have implemented this by making a linear interpolation of the loga¬ 
rithm of the data plotted against the logarithm of the momentum for both data 
sets. In this way the scaling of the momentum is accounted for by shifting the 
fine lattice data to the right by an amount A a as follows 

In D c (\n(qa c )) = lnD/(ln(ga c ) - A a ) + A z (13) 

Here A z is the amount by which the fine lattice data must be shifted up to 
provide the optimal overlap between the two data sets. The matching of the 
two data sets has been performed for values of A a separated by a step size 
of 0.001. A z is determined for each value of A a considered, and the optimal 
combination of shifts is identified by searching for the global minimum of the 
X 2 /dof. The ratios R a and R z are related to A a and A z by 

Ra = e” A “, Rz = R 2 a e~ Az . (14) 

Fig. I shows the data for both lattice spacings as a function of qa before 
shifting. In the leftmost plot of fig. [| we present the result of the matching 
using q as the momentum variable. The minimum value for y 2 /dof of about 1.7 
is obtained for R a ~ 0.815. This value for R a is considerably higher than the 
value of 0.716 ± 0.040 obtained from an analysis of the static quark potential 
inti From this discrepancy, as well as the relatively high value for x 2 /dof, we 
may conclude that the gluon propagator, taken as a function of q , does not 
exhibit scaling behaviour for the values of (3 considered here. 

The rightmost plot of fig. |] shows the result of the matching using q of 
(|[) as the momentum variable. We can see immediately that this gives much 
more satisfactory values both for y 2 /dof and for R a . The minimum value for 


4 





Figure 2: The dimensionless, unrenormalised gluon propagator as a function of ln(ga) for 
the two values of (3. The triangles denote the data for the small lattice at (3 = 6.0, while the 
circles denote the data for j3 = 6.2. 




Figure 3: x 2 P er degree of freedom as a function of the ratio of lattice spacings for matching 
the small and fine lattice data according to (|13j) , using q (left) and q (right) as the momentum 
variable. The dashed line indicates the ratio Rz of the renormalisation constants. 


X 2 /dof of 0.6 is obtained for R a = 0.745. Taking a confidence interval where 
X 2 /dof < Xmin + 1 givesJis an estimate of R a = O. 745 I 37 , which is fully 
compatible with the valueS of 0.716 ± 0.040. The corresponding estimate for 
the ratio of the renormalisation constants is Rz = 1.03811 1 - 

5 Model Fits 

We have demonstrated scaling in our lattice data over the entire range of q 2 
considered, and will now proceed with model fits. The following functional 
forms have been considered: 

D{q 2 ) = (\In ^±£) ~ dD (Gribov 3 ) (15) 
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D(q 2 ) = 
D(q 2 ) = 
D(q 2 ) = 
where 
D(q 2 ) = 
D{q 2 ) = Z 




Q 2 ( 1 l n ' 

+2A^V+A/ 4 ^ 2 111 AT 2 J 


(Stingfl)(16) 


z 

( 9 2)i + Q+M 2 


(Marenzoni ef all 


( g 2 +M 2 ( g 2 ))ln ’ a+ 4 ff 2(g2) 


(Cornwal. 




M(g 2 ) = Af {In 2^5^/In 

1 A 

9 2 +M^ v V 2 


4M‘ 

“a7 


■} 


- 6/11 


9 +M UV 


P) 


( 9 2)i+ q + ( m^ r )i+- ^ <J 2 +A/ 2 

D( g 2 ) = z (ie-f® 2 /"®)” + 




—cic 


(17) 

(18) 

(19) 

( 20 ) 
( 21 ) 


We have also 


considered0 


special cases of the three forms (Il9]) (Gil) , with 


M m = M uv or with specific values for the exponent a. Equations (1151) and (16) 
are modified in order to exhibit the asymptotic behaviour of ©• Models dT9[ ) 
and (||^) are constructed as generalisations of (0) with the correct dimension 
and asymptotic behaviour. 

All models were fitted to the large lattice data using the cylindrical cut. 
The lowest momentum value was excluded, as the volume dependence of this 
point could not be assessed. In order to balance the sensitivity of the fit be¬ 
tween the high- and low-momentum region, nearby data points within A (qa) < 
0.05 were averaged. 

X' 2 per degree of freedom and parameter values for fits to all these models 
are shown in table ||. It is clear that model (^o|) accounts for the data better 
than any of the other models. The best fit to this model is illustrated in fig. || 


6 Discussion and Outlook 


We have evaluated the gluon propagator on an asymmetric lattice with a large 
physical volume. By studying the anisotropies in the data, and comparing the 
data with those from a smaller lattice, we have been able to conclude that 
finite size effects are under control on the large lattice. 

A clear turnover in the behaviour of q 2 D{q 2 ) has been observed at q ~ 
lGeV, indicating that the gluon propagator diverges less rapidly than 1/q 2 in 
the infrared, and may be infrared finite or vanishing. 

The data are consistent with a functional form D(q 2 ) = Dir + D\jv, where 


-Dir 


1 1 

2 q 4 + M 4 ’ 


( 22 ) 
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Table 2: Parameter values for fits to models (|l5|-(|2l|). The values quoted are for fits to the 
entire set of data. The errors denote the uncertainties in the last digit (s) of the parameter 
values which result from varying the fitting range. The fitting ranges considered when 
evaluating the uncertainties are those with a minimum of 40 points included and with the 
minimum value for qa no larger than 1.03 (point number 40). 


M ~ 1 GeV, and -Duv is the asymptotic form given by ©• A more detailed 
analysis!^ of the asymptotic behaviour reveals that the one-loop formula (10) 
remains insufficient at q 1 2 3 4 5 = 50GeV 2 . 

Issues for future study include the effect of Gribov copies and of dynamical 
fermions. We also hope to use improved actions to perform realistic simula¬ 
tions at larger lattice spacings. This would enable us to evaluate the gluon 
propagator on larger physical volumes, giving access to lower momentum val¬ 
ues. 
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